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Radiative corrections to an atom are calculated near a half-space that has arbitrarily-shaped 
small depositions upon its surface. The method is based on calculation of the classical Green’s 
function of the macroscopic Maxwell equations near an arbitrarily perturbed half-space using a 
Born series expansion about the bare half-space Green’s function. The formalism of macroscopic 
quantum electrodynamics is used to carry this over into the quantum picture. The broad utility 
of the calculated Green’s function is demonstrated by using it to calculate two quantities — the 
spontaneous decay rate of an atom near a sharp surface feature, and the Gasimir-Polder potential 
of a finite grating deposited on a substrate. Qualitatively new behaviour is found in both cases, 
most notably in the latter where it is observed that the periodicity of the Gasimir-Polder potential 
persists even outside the immediate vicinity of the grating. 


I. INTRODUCTION 

Quantum fluctuations of the electromagnetic field are 
influenced by material boundaries, meaning that a wide 
variety of quantum electrodynamical vacuum effects have 
an environment-dependence which are often referred to 
as dispersion forces. Famous examples include the force 
between macroscopic objects known as the Casimir ef¬ 
fect [1], and the closely-related Casimir-Polder force [2] 
between an atom and a surface. Other examples include 
modified spontaneous decay rates [3, 4], magnetic mo¬ 
ments [5, 61, cyclotron frequencies [7, 8l and Zeeman 
splittings [9, 10]. 

There is contemporary interest in how dispersion forces 
are modified by the specifics of the surfaces involved. 
This can be by consideration of their optical properties 
[11], their thermal environment [12] or their geometries. 
An example of the latter is found in [13] where it is shown 
that nontrivial geometry-dependent vacuum effects can 
be studied by using a Bose-Einstein condensate above 
a corrugated surface. Dispersion-force calculations that 
go beyond simple planar geometries are usually compli¬ 
cated in the extreme due to the inherent non-additivity of 
dispersion forces (see, for example, [14, 15]). To remedy 
this, various simplifying approaches have been developed, 
one of the most prominent being the ‘proximity-force ap¬ 
proximation’ (PEA) [16] where one models complex ge¬ 
ometries as made up of an ensemble of flat, parallel sur¬ 
faces. It has been numerously shown (see, for example, 
[17-21]) that the PEA is uncontrolled and is often sig¬ 
nificantly in error. There are other approaches based on 
surface being ‘almost smooth’ [22, 23], but none are read¬ 
ily applicable to mechanically-etched surfaces with sharp 
edges like those discussed in [24], for example. 

Here we will use an alternative method based on the 
Born expansion of the Green’s function of the electro¬ 
magnetic wave equation, which will be used to calculate 
environment-mo diked decay rates and Casimir-Polder 
potentials near a selection of geometries. In contrast 
to the PEA, this approach preserves the rich geometry- 


dependence of dispersion forces, at the expense of requir¬ 
ing the system to consist of a small ‘geometric pertur¬ 
bation’ from an exactly solvable ‘background’ geometry. 
This approach has been used before in the calculation of 
Casimir-Polder potentials [25] and Casimir forces [26, 27]. 
One of the main differences between our work and [25- 
27] is that only homogenous backgrounds were considered 
there, while we consider a half-space as the background. 
The advantage of this is that the optical properties of the 
half-space can be specified completely freely — it is not 
part of the perturbation so its electromagnetic response 
does not need to satisfy any of the conditions that ensure 
convergence of the perturbation series. This allows one 
to make perturbative calculations for quantum electro¬ 
dynamical quantities near arbitrarily-shaped small depo¬ 
sitions onto the surface of the (non-perturbative) half¬ 
space, which is the goal of this paper. These kinds of 
geometries are relevant to very recent experiments on 
decay rates near pattered materials [28], and could also 
be applied to studies of surface roughness [29, 30]. 

II. THEORETICAL BACKGROUND 

We will use the noise-current approach [31, 32] to elec¬ 
tromagnetic (EM) field quantisation in and around di¬ 
electric media. This approach is necessitated by the 
fact that Maxwell’s equations in a dispersive, absorbing 
medium cannot be quantised simply by promoting the 
field observables to operators as this would cause a vio¬ 
lation of the fluctuation-dissipation theorem. To remedy 
this, one introduces a source current density operator j 
which corresponds to noise associated with loss in the 
medium and restores consistency with the fluctuation- 
dissipation theorem [31-34]. It is interesting to note that 
in its original form this theory did not rest on a rigorous 
canonical foundation, however this was recently reme¬ 
died in [35] . The advantage of the use of this source cur¬ 
rent representation is that it allows the quantised field to 
be obtained from the classical Green’s function for the 
electromagnetic field in a given configuration [31-34]. In 
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this framework, the frequency-domain quantised electric 
field that solves Maxwell’s equations in a medium with 
position and frequency-dependent permittivity e(r,co’) is 
given by the solution to the following wave equation [36] 

V X V X E(r, uj) — co’^e(r, co’)E(r, uj) = uj) , (1) 

with j being the operator-valued noise-current source dis¬ 
cussed above. This can be solved by the introduction of 
a Green’s function (variously called the dyadic Green’s 
function, or the Green’s tensor) [31, 32] which we will 
call W(r, r',cc;) [37]. It is defined as the solution to 

V X V X W(r, r',cj) 

— Co’^e(r,co’)W(r, r',cc;) = M(r — r'), (2) 


side the integrals, giving 

W(r,r',a;) = W(°)(r,r',w) 

Jw 

^ f d^si f d^S2 W^^^(r, si,u;) 

Jv iv *- 

X (si, S2, CJ)(S2, r', u;)l + .... 


( 4 ) 


In order to work out surface-modified quantities we will 
need the so-called ‘scattering’ part of the Green’s func¬ 
tion (4) — that is, the part which remains after the sub¬ 
traction of the Green’s function for a homogenous region. 
We write the scattering part of the whole Green’s func¬ 
tion W as G and the remaining homogenous part as H. 
This means the whole Green’s function can be rewritten 


where I is a 3 x 3 unit matrix. 

The Green’s function W defined by (2) uniquely de¬ 
termines the quantised field in a particular configuration, 
which ultimately means that knowledge of W allows one 
to calculate a wide variety of quantum electrodynami¬ 
cal quantities. However, exact calculation of the Green’s 
function W is only possible analytically for the very sim¬ 
plest choices of e(r,co’), so here we avoid this problem by 
using a perturbative technique. As shown in [25] one can 
write the unknown W in terms of some known ‘back¬ 
ground’ Green’s function (r, r', cj) as 


W(r,r',cj) = W^^^(r,r',cj) 

J (r, si,uj)Se{si,uj)W^^^ (si, r', uj) 

, Si, Ci;)^e(si, cj) 

X W*^^^(si,S2,cj)(5e(s2,cc;)W‘^^^(s2,r',cc;) -h , (3) 


where (5e(r,CL;) is the difference between the entire dielec¬ 
tric function and that of the background material at a 
particular point r. This type of perturbative expansion is 
known as the Born series and is the foundation of much of 
scattering theory — the spatial integrations over have 
a definite interpretation as scattering events [25-27, 34]. 

We can simplify the Born series (3) by specifying that 
the configurations we are interested in are always made 
up of an object described by some volume V that has an 
internally homogenous dielectric function e(r,ci;) = e(ci;), 
and sits in some (possibly inhomogenous) ‘background’ 
material with dielectric function e^^^(r,CL;). Under these 
assumptions we can restrict the integrals to being 
over the volume V, because outside this region the back¬ 
ground dielectric function at a particular point is equal to 
the entire dielectric function at that point, so (^e(r, cj) = 0 
there. The assumption of homogeneity within the volume 
V means we can also bring the dielectric functions out¬ 


W(r,r',w) = G(°)(r,r',w) +H(°)(r,r',w) 


+ UJ‘ 


'[(5e(w)] [ d^si 


G(°)(r,s,a;) + H(°)(r,s,w) 


G(°)(s,r',w) H-H(°)(s,r',w) 


+ ... 


( 5 ) 


In previous calculations [25-27] the background Green’s 
function W^^^ was taken to be that for a homogenous 
medium, so that its scattering part is by definition 
zero. This has the simplifying property that the parti¬ 
tioning of the Green’s function via the Born series coin¬ 
cides with the partitioning one makes when finding the 
scattering part, i.e. for homogenous H 

WHhom(r,r',u;) = (r, r', cj) 

uj^[6e{uj)] f d^SiH*^^^(r, s,co’)H^^^(s,r',cc;), (6) 

iv 

meaning that the scattering part is: 


GHhom(r,r',u;) = 

uj‘^[6e{uj)] f s,cc;)H^^^(s, . (7) 

Jv 

This means that all that is required for the calculation 
of environment-dependent quantities in a geometry re¬ 
garded as a perturbation to a homogenous medium is an 
integral over a homogenous Green’s function, which is rel¬ 
atively simple to do. However this is not usually a case 
of physical interest since in real experiments there will 
likely be an object nearby that does not obey the condi¬ 
tions for convergence that the Born expansion requires. 
To remedy this, we will study the simplest inhomogenous 
background, namely a half-space, and then add perturb¬ 
ing objects to thah as shown in Fig. I. 

In this work we will truncate the Born series at the 
single-scattering term, though it is straightforward to ex¬ 
tend the method to higher-order terms as is required to 
work out Gasimir forces [27] as opposed to decay rates 
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FIG. 1. General setup 


and Casimir-Polder potentials as is done here. We then 
have the whole Green’s function to order (5e(r,cc;) 

W(r,r',a;) = W“®(r,r',a;) 

W[(5e(w)] [ (fsW^^{r,s,Lo)W^^{s,r',uj) 

Jv 

= W“®(r,r',a;) + AW“®(r,r',w) , (8) 

with being the Green’s function for a half-space. 

The half-space Green’s function at frequency cj in a re¬ 
gion ^ > 0 in the presence of a non-magnetic material 
half-space filling the region ^ < 0 is conveniently written 
as [38]: 


W^®(r,r',w) = - 

a 


^ w ^ c-/ 

P.(r,r') 

k?\kz 


-r') 


(9) 


where k^ = y ry and ky are respectively the 

components of the position and wave vector parallel and 
perpendicular to the interface, and z is a unit vector per¬ 
pendicular to the interface. The symbol a indexes the 
two possible polarisations [TE (transverse-electric) and 
TM (transverse-magnetic)] of the Coulomb-gauge elec¬ 
tromagnetic field, and Vo- represents the following differ¬ 
ential operators 

^TE(r, r') = (V X z) (g) (V' X z) , 

^tm(i*, r') = ^ (V X V X z) O (V' x V' x z) . (10) 

Finally, the function F^{z^ z') is given by 

F^{z,z') = , ( 11 ) 

where is the greater of 2 : and z' ^ and Z< is the lesser 
of 2 : and z'] 

\ z ioT z > z' \ z' ioT z > z' 

^ yz' ioT z < z' ' ^ ^ 2 : for 2 : < 2 :' ’ 

( 12 ) 


and are the Fresnel coefficients for radiation propa¬ 
gating from a vacuum region into a medium of permit¬ 
tivity e{(jj) 


te _ kz k^ _ e{uj)kz k^ 

~ kzFkf ~ e{uj)kz + ki ’ 


(13) 


where k^ = e{uj)uj‘^ — /cy is the 2 :-component of the 

wave vector inside the medium. We can now use this 
statement of the half-space Green’s function to gener¬ 
ate the next-to-leading order term AW^^(r, r', uj) in the 
Born expansion ( 8 ), which will give the modified Green’s 
function for the EM field in the vicinity of a half-space 
with depositions. 


III. MODIFIED GREEN’S FUNCTION 


We now present the Green’s function modification 
AW^S(^j.^ for a half-space with a deposition. We 

will restrict ourselves to the region r, r' 7 ^ s throughout 
this work, meaning that we can ignore the S function 
terms in (9) when substituting it into ( 8 ). This means 
that we will not calculate any quantum electrodynamical 
quantities inside a deposition onto a half-space. Apart 
from complicating the method used here, calculation of 
such quantities would require the use of local-field cor¬ 
rected Green’s tensors [4, 39, 40] which are beyond the 
scope of this work. Under these assumptions, we note 
that AW(r, r',cc;) depends quadratically on so from 
the form of Eq. (11) one sees that that all contributions 
to AW as defined in Eq. ( 8 ) must be at most quadratic 
in the refiection coefficients, so we can write 


AW; 


^^{r,r',aj) = j J d^k\. 




where, for later convenience, we have defined the quantity 
P as 


P = -- 


Se{uj) 


exp 


{^kii - (i’ll -®ii) 


64:7T^uj^k‘^^k'^^kzF^ 

+ k\^-{s\\-r'\\) + k^{r^+s^) + k'^{r'^+s^) j. (15) 


The various Kij^ in (14) are matrix elements determined 
from Eqs. ( 8 ) and (9) by simple but tedious application of 
the differential operators ( 10 ) to the functions F^{z^z') 
defined in Eq. (11). The matrix elements differ depending 
on wether Vz is greater or less than Sz^ the subscript ^ 
distinguishes these two cases, as detailed in the full list 
of matrix elements found in Appendix A. 
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A. Simple demonstration: Decay rate near a sharp 
surface feature 


which we split into the free-space contribution Fq and a 
surface-modified part AF^^: 


We will begin with a surface-modified quantity that 
is relatively easy to calculate, namely the spontaneous 
decay rate F of an excited atom that is attributable to 
its interaction with the quantised electromagnetic field. 
It has been shown [41] that this rate can be expressed in 
terms of the Green’s function W as: 


r = 2 wAd • [ImW(rA,rA,WA)] -d* , (16) 

where d is the dipole moment of the transition. As an 
example we will calculate the decay rate Fq of an atom in 
vacuum, with no material objects present. We can choose 
the direction of the polarisation freely because in vacuum 
we have rotation invariance — we choose the polarisation 
to be aligned along the z direction so that d = dz. Then; 

To = 2wA|dpImW™‘'(rA,rA,WA), (17) 

where W^^*^ is the 2 : 2 : component of the Green’s function 
that solves (2) for e(r, cj) = 1. This vacuum Green’s func¬ 
tion is well-known (see [34] for a thorough review). It can 
be found, for example, from the half-space Green’s func¬ 
tion (9) reported here by taking all reflection coefficients 
to zero; 

W™‘=(r,r',a;) = W«®(r,r',w)|i^.=o ■ (18) 


The 2 ) 2 ) component of the vacuum Green’s function is 


W;r(r,r',a;A) 



i ^2 

_ \\_Jk^\z-z'\ 


(19) 


^HS 


= 2uj\d- [lin(Wvac(rA,rA,WA) 


+ G^^(rA,rA,WA) 


= + 2a;Ad • IinG^®(rA, ta, wa) • d* 

OTT 

= To + . (23) 


We will consider the two cases that the atom is polarized 
parallel and perpendicular to the surface with the same 
magnitude of dipole moment d = |d|, and separately find 
the two contributions AF^^ and AF^^ to the decay rates 

F|| = Fo + AF^^ and F^ = Fq + AF^^. Galculation 
of ATf is simplified by exploiting invariance in the xy 
plane to assume without loss of generality that the dipole 
in this case is aligned along the x direction. Therefore 
we need to calculate 


Arp = 2wi|dpImGP(rA,rA,WA), (24) 

Arp = 2wA|dpImGP(rA,rA,WA). (25) 

Using the half-space Green’s function (9) we find: 


pOO 

G^!{r,r',UA) = — / 

47r Jo ^ 


>.HS 


si 


2ikzz' 


^ikz{z-z') 


X [lu^Rte — • (26) 


where we have ignored the (real-valued) 6 function part 
of (9) in anticipation of taking the imaginary part as 
dictated by (16). Transforming to polar co-ordinates in 
the kx , ky plane and doing the trivial angular integral we 
have 


The integral can be carried out analytically, giving: 


W:r(r,r',n;A) 


1 f,iuj\z-z'\ 

- _ (1 

27rco’^ \z — 2 )'|^ 


— iuj\z — z'\) . ( 21 ) 


Taking the imaginary part of this followed by the coinci¬ 
dence limit z' ^ z^ we find upon substitution into (17) 


To = A\dA ( 22 ) 

which is a well-known result (see, for example, [42]), and 
will be used as a convenient unit in the following discus¬ 
sions. 

As another point of comparison we also present the 
results for the decay rate near a half-space 

r«® = 2uld ■ [lmW«®(rA,rA,:uA)] ■ d* , 


Substituting these into eqs. (24) and (25) and evaluating 
the integrals in the same way as for the free space case, 
one eventually finds the following results the decay rates 
near a perfect conductor {Rte ^ — 1, ^tm ^ 1) 


AF«s = 


id |2 


IGttz^ . 


(l — 4uj\z‘^) sin( 2 cjA^) 

— 2ci;a^ cos(2co’a2)) 


(27) 


ApHS ^ 


sin( 2 co’A^) “ 2 co’a 2 ) cos( 2 ci;a^)] , (28) 


in agreement with [43]. The 2 )-dependence of F^^ and 
F^^ is shown in Fig. 2. This shows the well-known prop¬ 
erty that an atom whose dipole moment is aligned per¬ 
pendicular to a perfectly reflecting surface has its decay 
rate enhanced by a factor of two in the small-distance 
limit. Similarly, an atom whose dipole moment is aligned 
parallel to such a surface has its decay rate completely 
suppressed as it approaches the boundary. Far away from 
the surface the free-space value is recovered in both cases 
as expected. 

We will now use our modified Green’s function (14) 
to produce new results for more complicated geometries, 
using the above known results as points of comparison. 
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FIG. 2. Decay rates near a simple half-space 



FIG. 3. Gubic deposition geometry. The cube and substrate 
can be made of different materials, the only restriction on 
their properties is that cube material must be weakly dielec¬ 
tric in order for the Born series to converge. We emphasise 
that the substrate can be made of any desired material since 
it does not take part in the perturbative approximation. 


The new geometry that we choose is a cube of side a 
and refractive index ec{uj) deposited on a half-space, as 
shown in Figure 3. This means the volume integral over 
s in Eq. (14) becomes 

p rCLl2 na 

/ d^s / dsx / dsx / dsz . (29) 

Jv J-a/2 J-a/2 Jo 

Part of the reason for choosing this shape in particu¬ 
lar is that, as mentioned in the introduction, the method 
presented here does not break down for geometries with 
sharp corners, in contrast to other approaches to radia¬ 
tive corrections near perturbed half-spaces which rely on 
the surface being smooth in some sense [22, 23]. As we 
will see later on, the approach used here can produce 
highly non-trivial results in the regions near sharp ob¬ 
jects. 

Taking the modified Green’s function (14) and 
transforming to polar co-ordinates {kx^ky} = 
{/c|| sin0,/c|| cos 0} (with similar for the primed co¬ 
ordinates) we find for the xx and 2 : 2 : components of 


the modified Green’s function in the limit where the 
substrate is perfectly reflecting: 

dk\\ dk\^ dtp dtp' 

k\\ ^11 kz COS (akz) sin {ak'^) — k'^ sin (akz) cos {ak'^) 

^ {k^ - kf){xk\\ - x'fc||)(?7fc|| - d'k'n) 


sm 


X -'n'k'w )+ry{xk\\ -Xk\\ )+/c^r^+/c^r; 

and 


“/2 (x^ll - X'fcy) 
“/2 {nk\\ -V^ii) 


(30) 




5ecM 

47r^ 


pOO pOO p2tt p2tt 

/ dk\\ / dk'\\ d(l) d(j)' 

Jo Jo Jo Jo 


kz sin (akz) cos{ak'^) — k'^ cos{akz) sm{ak'^) 

^ LO^k^k'^ {k^ - fcf )(xfc|| - X'k\{)ijlk\\ - 

X sin “/2 (x^ll - x'^) sin «/2 {dh “ V^) 

X k\^)+ry{xk\\-x'k\\)+k'yz+k^r^] ^ j^^l) 

where we have abbreviated 


X = cos 0, x' = cos 0', 


Tj = sm0. 


rj' = sin0'. 


(32) 


and immediately taken the parallel coincidence limit 
rj| ^ ry. The perturbative approximation holds as long 
as ec(co’A) — 1 < 1, meaning that we require ec(cJA) < 2. 
In practice this means that the absorption lines of the 
medium constituting the cube must be well-separated in 
frequency from the relevant atomic transition frequency 
uja- Here we will simply choose ec(co’A) to be such that 
the condition ec{oJA) < 2 holds. 

The quadruple integrals (30) and (31) are straight¬ 
forward to numerically evaluate in ready-made software 
such as Mathematica or Maple — no specialised numer¬ 
ical techniques are required. Their ease of evaluation 
arises because the angular integrals are over a finite range 
and the k\\ integrals are exponentially damped at infin¬ 
ity. We present a selection of results of this numerical 
study in Figs. 4 and 5. We note in particular that 
Fig. 5 shows the highly non trivial position-dependence 
of the decay rate — for example the decay rate can be 
enhanced or suppressed (relative to the value near a bare 
halfspace) depending on the precise position of the atom 
in the plane above the cube. 


IV. CASIMIR-POLDER POTENTIAL OF A 
FINITE GRATING 

A. Background and motivation 

We now turn our attention to a more complex but 
experimentally-relevant situation, namely the Casimir- 
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FIG. 4. Modified decay rates (solid lines) near a cube of di¬ 
electric constant ec(coA) = 1.8 deposited on a perfectly reflect¬ 
ing half-space, with the decay rates for the bare halfspace and 
free space (dashed lines) shown for comparison . The shaded 
area represents the depth of the cube added to the halfspace. 
We do not present results for the region interior to the cube 
because local-field effects [4, 39, 40] would come into play 
there, but these are beyond the scope of this work. 



X {CUaZ/c) 

FIG. 5. Normalised decay rate for an x-polarised atom at 
a distance 0.01a above a cube of side length a = 1 (in di¬ 
mensionless units ujAzIc) and refractive index eduj) = 1.8 
deposited on a perfectly reflecting half-space. The decay rate 
is expressed in units of the decay rate at the same distance 
above the ‘bare’ perfectly reflecting half-space (i.e. that with 
no cubic deposition). 


that for the decay rate in section III A. As we shall se, this 
is largely because the potential depends on a sum over 
all photon frequencies, rather than being determined by 
a specific transition frequency like the decay rate. An 
additional complication is that calculation of a CP po¬ 
tential involves subtraction of the contribution of the ho¬ 
mogenous part of the Green’s function at each particular 
point in order to extract a geometry-dependence. This is 
necessary because, unlike the decay rate, evaluation of a 
CP potential in free space (i.e., the Lamb shift) requires 
a completely different full field-theoretic approach. As 
detailed in the introduction, care must be taken with CP 
potentials in this Born-series approach because of the in¬ 
terplay between this subtraction of a homogenous part 
and the perturbative approximation. 

We will calculate the CP potential in vacuum near an 
N-grooved finite grating, like that shown in Fig. 6. This 
choice is motivated by the structures used ongoing ex¬ 
periments in atom optics and matter-wave interferometry 
such as [24] and [46]. There is section of existing litera- 



FIG. 6. Geometry of the finite grating considered here. Just 
as for the cubic deposition shown in fig. 3, the substrate can be 
made of any desired material, again the only restriction on the 
system is that the grating material must be weakly dielectric. 
While our method works for any number of grooves A, we will 
choose A = 5 here and throughout as shown in the figure. 

ture on CP forces near periodic gratings [17, 47], however 
these works take advantage of the Bloch theorem and so 
are only strictly applicable to infinite, precisely periodic 
gratings which are not necessarily good approximations 
to real experiments. In fact, as we will see later, non¬ 
trivial behaviour of the CP potential occurs outside the 
immediate vicinity of the grating, which of course cannot 
be seen if the grating is assumed to be infinite. 


Polder (CP) potential of an atom near a surface, as first 
described in [2] . The CP potential results from the mod¬ 
ification of the level structure of a polarizable atom by 
a surface-dependent quantised field — it is the surface- 
dependent version of the Lamb shift. The resultant force 
has been measured to high precision [44] and is of in¬ 
creasing importance in emerging quantum technologies 
[45] . The calculation is inherently more complicated than 


B. Expressions for Casimir-Polder potential 

The CP potential Ucp for an isotropically polarisable 
atom at position ta = (^a, Va^za) in a region with scat¬ 
tering electromagnetic Green’s function G may be writ¬ 
ten in terms of an integral over complex frequency ^ as 
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[34, 48] 

1 

Ucp{rA) = ^ j , (33) 

where a is the polarisability of a ground-to-excited 
atomic state transition of frequency uoij and dipole mo¬ 
ment dij and is given by 


/X 2 ujij dij 

a{uj) = - hm ^ 


3 uj‘f. — uj‘^ — iuje ’ 


(34) 


where 5 is a real infinitesimal [49] related to the line width 
of the atomic state [50]. The Green’s function is to be 
taken with both spatial arguments equal to the position 
ta of the atom, this is to be understood as a limiting 
value. Just as in the decay rate calculation in section 
III A, we will use a selection of standard results as points 
of comparison for later results. The first of these is the 
Casimir-polder potential at a distance from a perfectly 
conducting plane in the non-retarded regime, This regime 
is where the round-trip time for a photon to travel from 
the atom to the surface and back is much smaller than the 
timescale associated with the atomic frequency. In other 
words it is the small-distance approximation if the atomic 
transition frequency is assumed to be a fixed constant. 
The well-known result in this regime is [2]; 

77cpo(wA^A«l) = -i^. (35) 


The second quantity we will use as a comparison is the 
force that an atom in this potential experiences, namely 


Tcpo(^a^a < 1) 


-A-^CPo(^A^A < 1) 


\di,? 


(36) 


which a statement the well-known Casimir-Polder force 
of attraction between a polarizable atom and surface, in 
this case in the non-retarded regime and for a perfectly 
conducting material. 

Equation (8) tells us that the Green’s function that 
encodes the behaviour of the EM field near the grat¬ 
ing is given by the sum of two terms: which de¬ 

scribes the unperturbed half-space and AW^^(r, r', cj) 
which describes the correction resulting from deposition 
of the grating on its surface. The GP potential (33) re¬ 
quires the use of a scattering Green’s function. Since 
the whole Green’s function W is linear in its two con¬ 
tributions and AW^^(r, r', cj) it suffices to find 

the scattering parts and AG^^(r, r', cj) of these two 
contributions separately, which together give the scatter¬ 
ing Green’s function G. Also, the linearity of the GP 
potential in the scattering Green’s function means that 
we can find the contributions from the two scattering 
parts separately. We will confine ourselves to the most 
physically-relevant region z^z' > 0, meaning that the 


subtraction of a homogenous part is achieved by setting 
all reflection coefficients in the Green’s function to zero 
and subtracting the resulting quantity. Eor the unper¬ 
turbed part G^^ we have: 

G“®(r,r',w) = W*^®(r,r',w) - Wg®^o(r,r',w) , (37) 

where W^^^Q(r, r', cj) coincides with the Green’s func¬ 
tion of free space since the region near the grating as 
assumed to be vacuum. 

The scattering part of the Green’s function correction 
(14) is obtained by subtracting the portion that remains 
when setting all reflection coefficients to zero. Gonse- 
quently, isolating the scattering part AG^^(r, r', cj) of 
AW(r, r',co’) is trivial because of the way it is stated in 
Eq. (14) — all one needs to to is remove the term inde¬ 
pendent of reflection coefficients, giving: 


AGf 


= J^d^s j J 


+ KffRl) 


, (38) 


with matrix elements Kij listed in Appendix A. Now that 
we have the Green’s function correction (38) we can find 
the correction to the CP potential resulting from the de¬ 
position of the grating on the half-space from 


1 

MJcpivA] = ^ J c^CC^a(*C)Tr AG(rA,rA,iC) , 

(39) 


C. Grating results and discussion 

The volume integral describing the A-grooved grating 
shown in Eig. 6 is: 

n N—1 pxo-\-(2n-\-l)w rL/2 ph 

I d s y ^ ^ / dsx I dsy I ds^, •, 

dv ^_Q d xq-\-2tiw d — Z//2 do 

(40) 

with xo = —w{N — 3 / 4 ) if the grating is such that the 
centre of the base of the middle groove is at = 0. 

Eor simplicity the half-space will be taken as perfectly 
reflecting, and the grating as non-dispersive €g{uj) = eg 
with eg — 1 < 1. We choose A = 5, which corresponds 
to the grating shown in 6. The atom’s polarisability is 
taken to be isotropic. Using the volume elements (40) the 
integrals over s in (39) become elementary. This leaves 
integrals over /cy, /cy, ^, 6> and 0' which may be evaluated 
numerically. Just as for the decay rate calculation in sec¬ 
tion III A, the integration is significantly simplified by 
the fact that the integrals over 0 and 0' are both over the 
finite region and the remaining integrals are all 

exponentially damped. A selection of results are shown 
in Eigs. 7 and 8 The results for the CP potential di¬ 
rectly above the grating show qualitative agreement with 
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FIG. 7. Non-ret aided Casimir-Polder potential above the 
grating shown in Fig. 6. The potential is shown in units of the 
perfect reflector potential (35), which is the potential if the 
grating had not been deposited on the half-space. All length 
scales in the problem can be expressed in terms of a reference 
length which cancels out (in the non-retarded regime) when 
normalizing to the bare perfect reflector result. Thus, the 
units on both axes of the above plot are arbitrary — in other 
words the result remains valid whatever unit is assigned to the 
X and y axes as long as the non-retarded approximation holds. 
The parameters describing the grating are /i = l,ie = 1 and 
L = 5 in the same units as the x and y axes. The dielectric 
constant is tg = 1.8. Almost invisible in this figure is a sug¬ 
gestion of interesting behavior ‘outside’ the grating along the 
X axis, this is shown in detail in Fig. 8. 


the infinite grating considered in [17], where it was ob¬ 
served that the potential is reduced between the grooves 
and enhanced above them, as compared to the planar 
result. However our results are not directly quantita¬ 
tively comparable with [17] due to the choices of materi¬ 
als made there not being consistent with our perturbative 
expansion. Our results for the region ‘outside’ the grat¬ 
ing were of course not seen even qualitatively in [17] due 
to that work’s assumption of an infinite grating. Here we 
have relaxed this assumption, and found that the period¬ 
icity of the Casimir-Polder potential continues laterally 
past the end of the grating, which to our knowledge is a 
previously-unseen phenomenon. 


V. CONCLUSIONS 




X (arb.) 


FIG. 8. Upper plot: Potential Ucp near the grating shown 
in Fig. 6. As in Fig. 7 the potential is plotted in units of its 
value Ucpo [Eq. (35)] near a simple planar perfect reflector, 
and all the parameters are the same as for Fig. 7 which is 
reproduced as an inset and schematically on the horizontal 
axes. Lower plot: Lateral (x-directed) force Fcp near the 
grating shown in Fig. 6 for the same parameters used in Fig. 7. 
In the same spirit as the other plots, we normalise to what 
the force would have been if the grating were not present, but 
the lateral force without the grating would of course be zero. 
For this reason we use the (constant) perpendicular force (36) 
as a unit instead. 


of an atom near a cube deposited on a half-space, find¬ 
ing the rich position-dependence shown in Fig. 5. Fi¬ 
nally we presented the Casimir-Polder potential of a fi¬ 
nite grating deposited on a substrate and demonstrated 
the previously-unseen quality that the lateral periodic¬ 
ity of the potential can continue beyond the grating, as 
shown in Fig. 8. The Green’s function (14) can be used 
to calculate quantum electrodynamical quantities near a 
half-space with any small deposition on it, so the work 
presented here should have applications in a variety of 
ongoing and planned experiments [24, 45, 46]. 


In this paper we have considered some aspects of 
quantum electrodynamics near a surface with arbitrar¬ 
ily shaped features deposited on it. The main general 
result is the Green’s function (14) for the perturbed half¬ 
space, which was calculated perturbatively using a Born- 
series expansion. We then investigated the decay rate 
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Appendix A: Green’s function matrix elements 

There are five terms in Eq. (14), each of which is a one 
of two 3x3 matrices (one for each of choice of ^ 
giving a total of 90 matrix elements that we, in principle, 
need to calculate. However there are various constraints 
that reduce this number significantly. Firstly, some ma¬ 
trix elements are not independent due to the xy symme¬ 
try of the half-space. In particular; 

~ ^ ky), 

^zy^ ~ ^zx^{kx ^ ky)^ 

K;x^ = Ky^^ (Al) 

where r = (TE, TM, TETE, TMTM, TETM). This 
restriction reduces the number of required matrix ele¬ 
ments by 

(4 constraints ) x (5 different r) x (2 for Vz ^ Sz) = ^0, 

leaving a total of 50. This number can be further reduced 
by noting that the definition of TE modes is that they 
have no electric field in the z direction, which ultimately 
means that any matrix element for TE polarisation where 
at least one index is z is in fact identically zero 


F^{z^z') given by Eq. (11). For Vz > Sz the matrix 
elements representing coefficients of terms linear in the 
reflection coefficients are: 

(klkl + ky) + primed, 

^Iy> = kxk'^kykyUj‘^ ^ + primed, 

[klkl + ky) + primed, 
A"™ = e~‘^^^^^^k‘^\k'^kykykxk^ + primed, 

A"Ji> = ^~‘^^^''^''kxk'^kzk'yy — primed, 

— ~^~‘^^^^^^kxk'^kzk'yy — primed, 

K™ = + primed (A4) 

where ‘primed’ is a shorthand for the quantity that pre¬ 
cedes it with k ^ k' and r ^ r', with the latter re¬ 
placement only being relevant for Vz < Sz as we shall 
see. Continuing, the Vz > Sz matrix elements represent¬ 
ing coefficients of terms quadratic in particular reflection 
coefficients are 


= k‘^,k[Juo^, 


TETE 

xx> '"y'"y' 

= k',kMyW\ 


K, 

K _ 

-‘^^xy> f^x'^x'^y'^y^ 

^TMTM _ 7 7/7 If 7^27/2 

^xz> '^x'^x^zf^z^yW 5 

^TMTM _ 7 4 7 /4 
^ZZ> — 5 


^TMTM _ 7 2^/2 7 2 7 /2 
^^XX> ~ '^x'^x '^z'^z 5 

^TMTM _ jL U y iu2 7^'2 
^^xy> ~ '^x'^x'^y'^y'^z'^z 5 

^TMTM _ 7 ifi 7 / 7^27/2 

^ZX> '^X'^x^z’^z^W^W 5 

(A5) 


t^TE _ ^TE _ t^TE _ Q 
^xz^ - ^zx^ - ^zz^ - 

t^TETE _ ^TETE _ ^TETE _ n 
^xz^ - ^zx^ - ^zz^ - ^ 5 


and similarly 


t^TETM _ t^TETM _ t^TETM _ n 
^xz$ — ^zx^ — ^zz$ — ^ 


(A2) 

(A3) 


which together reduce the required number by 18 leaving 
90 — 40 — 18 = 32 matrix elements to calculate, which can 
be partitioned into two groups of 16, where each group 
corresponds to one choice of rz ^ We now simply 
list these matrix elements, which are obtained by appli¬ 
cation of the differential operators (10) to the functions 


and Anally the coefficients of the terms that mix TE and 
TM reflection coefficients 

= {kXkl + k'ykf), 

= kykyk’y {kl + kf) . (A6) 

For Vz < Sz the entire set of 16 coefficients can be ob¬ 
tained from Eqs. (A4)-(A6) by taking Sz Vz (before 
adding the ‘primed’ parts), so that for example: 

K™ = + e-2"=-4)qfc'4^ (A7) 

We have now completely specified all terms in (14). 
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